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Abstract 

We apply the method of Dalgarno and Lewis to scattering states 
and discuss the choice of the unperturbed model in order to have a 
convergent perturbation series for the phase shift. 

1 Introduction 

There has been some interest in the application of on-shell perturbation the- 
ory to simple scattering problems in one dimension [1-3]. Those approaches 
were mainly based on logarithmic perturbation theory [1] and the Dalgarno- 
Lewis perturbation theory [2,3]. In order to bypass the problem posed by 
the nodes of real wavefunctions they used complex wavefunctions (Jost func- 
tions). Although the method of Dalgarno and Lewis was proposed as a device 
to calculate sum over states [4] it proves to be an efficient approach for on- 
shell perturbation theory [5,6]. This method has recently been applied to 



solve the differential equation for the density distribution in scattering prob- 
lems [7]. 

The purpose of this paper is to discuss the application of the Dalgarno and 
Lewis perturbation theory to potential scattering using real wavef unctions. 
We also consider the choice of the unperturbed or reference model in order to 
have convergent perturbation series for the phase shift. In Sec. |2]we introduce 
the model. In Sec. |3]we propose an alternative application of the method of 
Dalgarno and Steward and a choice of the unperturbed model potential. In 
Sec. |Uwe apply the main equations to an illustrative example, and in Sec. |H1 
we summarize the main results of the paper and draw some conclusions. 

2 The model 

For simplicity we consider the one-dimensional model 

h 2 

- ^—^"(r) + V{r)ij){r) = Eip(r) (1) 

and choose the boundary condition ip(0) = in order to mimic the s-states 
of a central-field model. The prime denotes differentiation with respect to 
the coordinate. The potential-energy function V(r) is negative for r < L 
and zero for r > L. In order to derive a dimensionless differential equation 
we define the new independent variable x = r/L and rewrite Eq. as 

y"(x) = [v{x) - e]y(x) (2) 

where 

2 772 2/77, 

v ( x ) = jf v ( Lx )i e = jf E i y( x ) = i>(Lx). (3) 

The boundary conditions are 

y(0) = 0, y(x) = ~ sin(A;x + 8), x > 1 (4) 

K 

where 5 is the phase shift and k = y/e. In order to avoid ambiguities regarding 
to multiples of n we choose 8(k) to satisfiy 6(oo) = 0. 

In order to apply perturbation theory we choose a closely related problem 

y'o( x ) = NO) - e]yo(x) (5) 
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with similar boundary conditions 

yo(O) = 0, y (x) = i sin(kx + 6 ), x > 1. (6) 

As it is custommary in most applications of perturbation theory we write 
v(x) = vo(x) + Xvi(x), where the perturbation parameter A is set equal to 
unity at the end of the calculation, and expand 

00 

y(x) = ^2yj[x)X 3 , 
3=0 
00 

8 = J^A'. (7) 

j=Q 

3 The method of Dalgarno and Lewis 

For simplicity we write Q{x) = v(x) — e, and Qq{x) = v (x) — e. The function 
F(x) = y(x)/yo(x) satisfies the differential equation 

(y {xfF\x))' = AQ(x)F(x)y (x) 2 (8) 

where AQ = Q — Qq. Since the nodes of y(x) and yo(x) do not coincide one 
expects F{x) to have poles at the nodes of yo(x). However, F(x) does not 
appear in the resulting expression for y(x) (see below) and, consequently, that 
problem does not arise in practical applications of the method of Dalgarno 
and Lewis. 
Notice that 

y (x) 2 F'(x) = y'(x)y (x) - y(x)y' (x) = W(y,y )(x) (9) 

is the Wronskian of the perturbed and unperturbed solutions. On integrating 
Eq. (jHJ) we have 

W{y,y )(x) = [ X (AQyy )(x')dx'. (10) 
Jo 

At x = 1 we substitute the asymptotic forms of y(x) and yo(x) into the 
Wronskian and obtain an expression for the phase shift 

sin(A5) = -k I (AQyy )(x)dx. (11) 
Jo 

3 



If we integrate Eq. (|1U|) between x = 1 and x and multiply the result by 
yo(x) we obtain an expression for y(x): 

y(x) = Cy (x) + y (x) / -— ^ / (AQyy ) (x") dx" (12) 

Ji y {x'y Jo 

where C is an integration constant. Notice that this expression is free from 
poles because the function F(x) does not appear explicitly in it. 

When x > 1 then x' > 1 and we substitute Eq. (jll)) for the second 
integral in Eq. ()12|) . Moreover, in order to have the correct asymptotic 
expression for y(x) the integration constant should be 

C = sin(A<5) cot(fc + 6 ) + cos(A5) (13) 

that we also expand in a power series: C = 1 + C\ \ + C2A 2 + . . . 
Finally, it follows from equations (|11|1. ()12|) . and ()13|) that 

/•a; ^2;' rx 1 

Vj{x) = Cjy (x) + y (x) / («i^_iy ) (z") dz" ( 14 ) 

<5i = —k / wi(x)yo(^) 2 ^ 
Jo 

<5 2 = -k v 1 (x)y (x)y 1 (x) dx 
Jo 

5 3 = -}- - k / wi(x)?/o(a;)l/2(a;) <ix (15) 
o Jo 

and 

d = 5i cot (A; + 5 ) 

C 2 = 5 2 cot(A; + 5o)-y 

C 3 = (s 3 -^jcat(k + 5 )- 5 1 5 2 . (16) 

At the jth perturbation step we first solve for yj(x) and Sj and then calculate 
Cj. Those expressions are sufficient for present discussion, one can easily 
derive as many as necessary from the equations above. 

The number of necessary perturbation corrections depends on the con- 
vergence rate of the perturbation series. For that reason it is important to 
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choose a convenient unperturbed or reference potential v (x). It has been 
argued that in order to have an adequate rate of convercence the unper- 
turbed and perturbed potentials should support the same number of bound 
states [1]. This conclusion is based on Levinson's theorem than in our case 
takes the form 5(k = 0) = Ntt, where N is the number of bound states [8]. 
Here we investigate the application of the Bargmann-Schwinger upper limit 
to the number of bound states [9, 10] and require that at least 

/ vo(x) dx — / v(x)dx. (17) 
Jo Jo 

Although this condition does not completely guarantee that both potentials 
have the same number of bound states, it is a reasonable choice that will 
prove sound in the example below. 

The simplest exactly solvable unperturbed model is one with a constant 
potential 

/ \ f — if < x < 1 , 10 s 
^ = ( if x > 1 ^ 

in which case ^ 

vq — — v(x) dx. 
Jo 

The starting-point of the calculation are the unperturbed solutions 

. , sin(/c + 5n) 

= u • r^- MKx), 0<x< 1, 
fcsm(A) 

K = ^k 2 + v , 

fktan(K)\ , 
o = arctan I — I — k + nir (19) 

where we choose n so that S (k — > oo) = as stated above. 



4 Example 

In what follows we compare the perturbation results with accurate numerical 
phase shifts calculated by a straightforward power-series method. To this 
end we assume that the Taylor series of the potential around x — 0, v(x) = 
wo + wix + . . . converges at x — 1, and expand the solution in the same way: 
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y{k,x) = c (k) + ci(k)x + c 2 (A;)a; 2 + . . ., where c = c 2 = and Cj(k,Ci 
Cj(k, l)c\, j = 3, 4, . . .. Finally, we obtain the phase shift as 

y(k,i 



arctan 



— k + H7l. 



(20) 



y(M), 

As an illustrative example, we choose 

u(ac) = Ac(x - 1), A > 0, < x < 1 (21) 

that can be solved in terms of confluent hypergeometric functions [11]. How- 
ever, here we resort to the power series approach because it is suitable for 
more general problems and converges fast in most cases. The power-series 
approach also enables us to obtain the bound-state energies from the condi- 
ton 

y'{k, 1) + ky{k, 1) = 0, k = V=~e (22) 

We choose a well of depth vq as unperturbed model and the simple con- 
dition (|T7|) gives us v = A/ 6. We could obtain analytical expressions for 
5i(A, vq, So), S2(A,v ,6 ) and S 3 (A,v ,6 ) but we only show the first pertur- 
bation correction because the other ones are too long. 

sin(5 + k) 2 



6i 



12k(k 2 + w ) 3 / 2 sin Uk 2 + v 



2v (k 2 + v 



sin 



k 2 + V ) COS A/ k 2 + Vq 



+ ^k 2 + v 



3a sin (yj k 2 + v ^j — a (k 2 + v + 3 



+6w ( k2 + v o)]} 



(23) 



In order to test the accuracy of perturbation theory we choose parabolic 
wells with A = 6 with no bound states, and A = 18 with one bound state. 
Figures 1 and 2 show the logarithmic error ei og = log\(8ps — 8pt)/Sps\, 
where PS and PT denote power series and perturbation theory, respec- 
tively. We clearly appreciate that the perturbation series converges for both 
parabolic-well strengths, and for all energies. For some k values the per- 
turbation series of order j — 1 may be more accurate than the one of order 
j. We do not know if this behaviour is fortuitous or if there is a mathe- 
matical or physical reason behind. Table ^ shows the average relative error 
e av = jjJ2fU\[Sps(kj) - 5 PT (k j )]/5p S (k j )\ for some A-values. We clearly 
appreciate that the perturbation series converges smoothly. 
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5 More general problems 

One can easily apply the method of Dalgarno and Steward to more general 
problems of the form 

y"(x)=Q(y,x) (24) 

where Q(y,x) denotes a differential operation on y(x). We choose a conve- 
nient unperturbed problem 

y'o(x)=Q (y ,x) (25) 
and proceed as in Sec. El We obtain the following general equation 

r x dx' 



y(x) 



C 2 yo(x) + Ciy (x] 
dx' 



+yo{x) 



yo(x 



A2 



p yo(x') 2 
{Qyo - Qoy) {x") dx 1. 



(26) 



where we can choose the arbitrary integration limits a and /3 and constants 
Ci, and C2, conveniently according to the problem. 
Notice that if Q (y ,x) = Qo(x)yo(x) then 



u(x) = y {x) 



x dx' 



p yo{x 



A2 



(27) 



is a solution of u"{x) = Qo(x)u(x) and W(u, yo) = 1. That is to say, yo(x) 
and u(x) are two independent solutions of the unperturbed model. 

It is not difficult to verify that the method of Dalgarno and Lewis applies 
to velocity-dependent problems [13]. If we choose 



Q(y,x) 



[v (x) + Xv^x) - e)y(x) + Xp(x)y"(x) 
y(x) 



-X 



y [x) 



p'(x) 



then we have the model treated by Jaghoub [13]. Therefore, if 

Qo(y ,x) = [v Q (x) - e ]y {x) 

we have 

Qyo - Qoy = 



(28) 



(29) 



-A 



Aeyy + Xvxyyo + Xpy y" 

y 



y - 



X 



p ya- 



rn 



10" 



-6 



o 



4 



k 



Figure 1: The error e 



<>PS 



x 100 for A = 6. 



If we expand y(x) as in Eq. (JZj) and Ae = eiA + e 2 A 2 + . . . then we obtain 
Jaghoub's perturbation equations [13] provided that we choose the arbitrary- 
constants conveniently. In other words: Jaghoub's procedure is merely the 
method of Dalgarno and Lewis for bound states developed in a different way. 

By means of the equations displayed above, the reader can easily verify 
that present implementation of the method of Dalgarno and Lewis is also 
suitable for the treatment of scattering problems with velocity-dependent 
interactions [12]. 

6 Conclusions 

We have shown that the method of Dalgarno and Lewis is suitable for the 
calculation of phase shifts by means of perturbation theory. One can use real 
eigenfunctions because their nodes do not cause the occurrence of poles into 
the solutions. 

The constant potential is a suitable unperturbed or reference model for 
the application of perturbation theory to some problems. A simple prescrip- 
tion for the well-depth leads to reasonable convergence rate and enables one 
to obtain accurate results with perturbation series of low order. 
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Table 1: Average relative error for the parabolic well V(x) = Ax(x — 1) 



.4 


Zeroth order 


First order 


Second order 


Third order 


6 


0.1060 


0.00736 


0.00024 


0.000018 


12 


0.0917 


0.02250 


0.00106 


0.000672 


18 


0.0906 


0.01260 


0.00127 


0.000365 


24 


0.0814 


0.01109 


0.00200 


0.000327 



We have shown that the method of Dalgarno and Lewis also applies to 
bound-state and scattering problems with velocity-dependent interactions 
[12,13]. 
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